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We propose an original adaptive wavefront holographic setup based on the 
photorefractive effect (PR), to make real-time measurements of acousto-optic 
signals in thick scattering media, with a high flux collection at high rates 
for breast tumor detection. We describe here our present state of art and 
understanding on the problem of breast imaging with PR detection of the 
acousto-optic signal. 
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1. Introduction 

The field of acousto-optic imaging has been strongly 
stimulated by the deep and complete paper of W. Leutz 
and G. Maret 1 . In this paper the authors give a very 
clear view of the tricky interactions between light and 
sound in random media; this is why this work has 
stimulated a new active field now more "biomedical 
imaging" oriented and new detection schemes. 

The present paper is at the frontier of two physical 
domains that are 

1. detection of weak light signal by using photorefrac- 
tive crystals, 

2. breast cancer imaging by detection of the ultra- 
sonic modulation of the light scattered through the 
breast. 

Here, our purpose is to make a brief review of these two 
domains, and to describe the photorefractive detection of 
the scattered light modulated component in a pedagog- 
ical manner. By the way, we will describe our present 
state of art and understanding on the problem. 

In this paper, we will first describe the basics principle 
of ultrasonic modulation of light imaging. We will, in 



particular, introduce the concept of "ultrasonic tagged 
photons", which represents the weak signal to be de- 
tected. We will then describe how photorefractive adap- 
tive holography can be used to detect the tagged pho- 
tons. One must notice that all the groups working on the 
subject, exempt to us, do not detect the "tagged" pho- 
tons, but "untagged" ones. The tagged photon signal is 
measured indirectly, since the total number of scattered 
photons (tagged + untagged) does not depend on the ul- 
trasound. We will describe our technique and present our 
experimental results. In all these descriptions, we must 
not forget one difficulty that results from the decorrela- 
tion of the light that travels through breast organ. This 
effect known as speckle decorrelation is both due to the 
brownian motion of the scatterers, and to the breast in- 
ner motions (blood flow ...). In a typical in vivo situa- 
tion, with 4cm breast thickness, the "speckle decorrela- 
tion time" is in the 0.1 to 1ms range. It is thus necessary 
to match the so called "photorefractive response time" 
with the "speckle decorrelation time". This effect, which 
is huge in breast, is not present in most of the ultrasonic 
modulation test experiments, which are performed with 
breast phantoms like dead tissues or diffusing gels. Since 
the decorrelation affects considerably the detection sen- 
sitivity, it is quite difficult to evaluate the figure of merit 
of the different techniques that are proposed to perform 
breast imaging. We will see that our setup, which is 
able to detect both the tagged and untagged signal, is 
also able to measure the photorefractive time in situ, i.e. 
with the same setup, same laser powers and same sample 
geometry than for breast imaging experiments. To our 
knowledge nobody is presently able to perform ultrasonic 
modulation imaging though 4cm of breast tissues in vivo. 
Experiments are under progress and we hope to be able 
to reach this aim in a near future. 
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2. Acousto-optic imaging 



The combination of light and ultrasound to measure lo- 
cal optical properties through thick and highly scattering 
media is a tantalizing approach for in vivo imaging. It 
is an alternative solution to pure optical techniques for 
breast cancer detection. The use of light is motivated by 
its relative low absorption in the so called "optical ther- 
apeutic window" (700nm to lOOOnm), and by the exis- 
tence of optical contrasts between healthy and tumorous 
areas in this region of the spectrum. 

Light is highly scattered within biological tissues, mak- 
ing direct optical study of thick sample very difficult to 
perform. Light scattering is characterized by two length 
parameters, e.g the scattering length l S: and the light 
transport mean free path I*. The scattering length l s 
characterizes the memory of optical phase, and corre- 
sponds to the average distance that separates two scat- 
tering events. The light transport mean free path I* char- 
acterizes the memory of the light propagation direction. 
In tissues, l s is typically 50 to 100 /im, while /J is 10 x 
larger (0.5 to 1mm). Absorption of light is characterized 
by the absorption length Z a , which is in the 1 cm to 10 
cm range. Absorption strongly depends on the nature of 
the tissue (optical contrast). 

Because of scattering, direct imaging cannot be per- 
formed through more than a few millimeter thick sam- 
ples. Contrarily to light, ultrasound (US) beams are 
ballistic in biological tissues. US gives thus access to 
millimeter range spatial resolution in thick sample (up 
to 4cm) yielding the development of the acousto-optic 
imaging that combines optics and ultrasound^. 

A. Principle: the tagged photons 

Acousto-optic imaging is a hybrid technique, which com- 
bines, thanks to the acousto-optic effect, ultrasound and 
light. US are applied in the region of interest, within 
the thick scattering sample (see FigJT] a). They make 
the scatterers vibrate. A CW laser (frequency uol) illu- 
minates the sample. The vibration of the scatterers at 
the acoustic US frequency uoa (2 MHz typically) modu- 
lates the phase of the photons that are scattered by the 
sample. This is the so-called acousto-optic effect. 

The light exiting the sample contains thus different fre- 
quency components (see FigJT] b). The main component 
(the carrier) is centered at the laser frequency uol- It is 
related to the diffused photons, that do not interact with 
the US. The sideband components are shifted by the US 
frequency uol ± wa- The sideband photons, which result 
from the interaction between light and US, are called 
"tagged photons" (i.e photons tagged by the US). 

The weight of the tagged photons components depends 
on the optical absorption in the region of interest, where 
the US beam is focused. Acousto-optics imaging stands 
in detecting selectively the tagged photons. An image 
of the sample optical absorption can be then built-up 
in scanning the US over the sample. Note that one of 
the difficulty in living tissues results from the motion 
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Fig. 1. Principle of acousto-optic imaging, a) Motion of 
the scatterer at uja- b) Spectrum of the diffused light: 
carrier (ujl), an d sideband i.e. tagged photons (cjl±uja)- 



of scatterers {e.g blood flow) which broaden the carrier 
and sideband lines (see Fig{T]b). In vivo measurements 
through 4 cm breast tissues yield a broadening of 1.5kHz 
(Full Width at Half Maximum: FWHM^. 

B. State of the art for the detection of the tagged pho- 
tons 

Many techniques have been proposed to detect the tagged 
photons. Marks et al.^ investigated modulation of light 
in homogeneous scattering media with pulsed ultrasound. 
Wang et al.— »£ performed ultrasound modulated optical 
tomography in scattering media. Lev et al. studied scat- 
tering media in the reflection configuration 8 . Wang and 
Shen 9 developed a frequency chirp technique to obtain 
scalable imaging resolution along the ultrasonic axis by 
use of a one-dimensional (ID) Fourier transform. Lev et 
al. use a set of optical fibers coupled to a single photo- 
detector— that allows to work with samples, which 
decorr elate in time. Leveque et a h n i 12 i 13 performed par- 
allel detection of multiple speckles on a video camera 
and demonstrated improvement of the detection signal- 
to-noise ratio of ID images of biological tissues. The 
parallel detection has been still improved by Gross et al., 
which performs holographic detection reaching the shot 
noise sensitivity limit 14 , and by Atlan et al., which get 
resolution on the US propagation axis by using an holo- 
graphic pulsed technique 1 ^. 

All these methods exhibit two main limitations. First, 
the optical etendue (defined as the product of the detec- 
tor area by the detector acceptance solid angle) of the 
detection system is not optimum, since it is much lower 
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Fig. 2. Principle of the PR detection of the acousto-optic 
signal. PZT: ultrasound transductor, Ej: illumination 
optical field; Ep,Eprp tagged or untagged field scattered 
by the sample; PR crystal: photorefractive crystal; L: 
lens that collect the scattered field into the crystal and 
photodiode; Er: PR crystal reference (or pump) field. 



than the etendue of the tagged photons source. This 
etendue is the area of the sample (several cm 2 ) x the 
emitting solid angle (which is about 2tt since the light is 
diffused by the sample in all direction) . With a mono de- 
tector (photodiode) 5,8,10 the detection etendue is about 
A 2 . With a multi detector like a CCD camera±±>±^>i^ 
the etendue is 7VA 2 , where N is the is CCD number of 
pixel (N ~ 10 6 ). Even with a camera, the etendue of 
detection is about x 1000 lower than the etendue of the 
emission. 

The second problem occurs within living sample: the 
scatterers move, yielding in the frequency space a broad- 
ening of the tagged photons spectrum, as shown on 
FigHfc 4,5 . This effect corresponds, in the time space, to 
a decorrelation of the tagged photons speckle pattern. 
Since all the methods described above perform coherent 
detection, the bandwidth of detection is limited by the 
detector bandwidth. With camera (there is no problem 
of bandwidth with photodiode, but the etendue is much 
lower), the bandwidth is roughly equal to the camera 
image frequency ojccd, which is in general much lower 
(^ccd ~ 10... 100 Hz) than the tissue broadening (3 
kHz). It is still possible to work with fast camera (kHz), 
but in that case i) the camera quantum efficiency is lower 
(CMOS), and ii) the number of pixel N is limited, be- 
cause N x uo cc d is the flux of information to transfer to 
the computer, and this flux is limited (< 10 6 ...10 T s -1 ). 

3. The photorefractive (PR) detection of the 
acousto optic signal 

More recently has appeared a new tagged photons detec- 
tion technique that is based on the photorefractive effect 
(PR) and that is illustrated by FigEl The light that is 
scattered by the sample (Ep or Eut for the tagged or 
untagged field) is detected by a photorefractive detec- 
tor (PR crystal + photodiode PD) that is pumped by a 
reference field Er. 

Since the crystal and the photodiode might be quite 





Fig. 3. Principle of holography using the photorefractive 
effect. 
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Fig. 4. Method for tuning the frequency ujr of the ref- 
erence beam: our = uop + uja for example. AOM1 and 
AOM2 are acousto optic modulators. 



large (up to 1cm 2 ) and since the light is collected by 
a large Numerical Aperture (N.A. ~ 1) collecting lens, 
the photorefractive detection benefits of a high etendue, 
about 100 x larger than in a typical camera with N ~ 10 6 
pixels. We will see that the detection bandwidth is the 
inverse of the "photorefractive time" Tpr. We get for 
example, for a lW/cm 2 pump beam, 1/Tpr ~ 1 kHzi£*i£. 
This bandwidth, which is about 100 x larger than for 
a N ~ 10 6 CCD camera x , is within the range of the 
linewidth ( ~ 3 kHz) of the light scattered in vivo by a 
breast organ 4 . 

A. The volume hologram 

Photorefractive effect arises in materials that present 
both electrooptic effect and photoconductivity, which 
combination allows to transform a non uniform illumi- 
nation of the material into a spatial variation of the re- 
fractive index 18 . When illuminated by the interference 
pattern between an object and a reference beam, the ma- 
terial records a hologram, i. e the amplitude and phase of 
the object beam. This hologram is dynamic meaning 
that it can follow the interference pattern fluctuations 
slower than the response time Tpr of the material, also 
meaning that only slowly moving hologram are recorded. 

The first effect is the recording of the signal beam in- 
formation (phase and amplitude of the signal field Ep) 



1 we implicitly exclude here fast CMOS camera because of poor 
quantum efficiency and noise, and because of finite bandwidth 
for the data transfert form camera to computer 
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within the PR crystal under the form of local changes of 
the crystal refractive index n yielding a volume hologram 
(FigEK). Since the recording takes a finite time Tpr (0.1 
to 10ms in our device), the PR effect selects, within the 
signal beam, the components whose frequency is close or 
equal to the reference beam frequency ojr. 

Here, the large, almost flat field, reference beam (field 
Er) interferes with signal field Et onto the PR crystal. 
In order to select the tagged or untagged photons, the 
frequency ujr of the reference beam is made equal to the 
tagged or untagged photons frequency: ujr = ojl or 
ur = oo l respectively. To adjust ojr one can use for ex- 
ample two acousto-optic modulators AOM1 and AOM2 
(Bragg cells) as shown on FigHl With such a choice, the 
interference pattern EtE r of the reference beam with 
the selected photons beam varies slowly in time. The se- 
lected beam information can thus be grooved within the 
PR crystal volume hologram. 

In the case of a perfect monochromatic signal beam, 
the local variation of the hologram refractive index 5n 
is simply proportional to the modulation depth of the 
interference pattern (E T E* R ) / (\E T \ 2 + \E R \ 2 ). If a time 
modulation is added on the signal (e.g. amplitude or 
phase modulation), we have to take into account the fi- 
nite time Tpr needed to groove the hologram, and we 
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where ( ) TpR is the average over the grooving time Tpr, 
average which is defined by: 
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B. The diffracted beam Ed 



The second effect is illustrated by FiglSb. The reference 
beam (Er) is diffracted by the volume hologram yielding 
a diffracted beam (Ed)- The diffracted field Ed is simply 
proportional to the hologram refractive index changes Sn 
and to the reference beam field Er. We get thus: 



E D oc 



(EtE 



R/T PF 



\E 7 



\E h 



,E R 



(3) 



In typical application the reference beam intensity is 
much larger than the signal beam one, and except of the 
average over Tpr, Er and E R simplifies in EqJH yielding 
E D oc E t i.e. 



E r 



(4) 



where n = 0.1.. 0.5 is a numerical factor which mainly 
depends on the crystal. 

EqHlis valid, when the decorrelation of the signal field 
Et can be neglected during the grooving time Tpr, i.e. 
when 



SojTpr <C 1 



(5) 
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Fig. 5. a,b) Detection of the tagged photons, when the 
phase of US beam is zero (a) and tt (b). c) Phase of the 
US beam. 



where Sou is the frequency width of the signal beam 
(Aoo ~ 3 kHz for the breast). 

We have to notice that an increase of the reference 
beam intensity \Er\ 2 does not change 77, but reduces the 
grooving time Tpr, since Tpr oc 1/\Er\ 2 . The main 
advantage of increasing the reference beam power \Er\ 2 
is thus to reduce Tpr enough to neglect the signal field 
decorrelation. Condition of Eqj5]is then fulfilled, and the 
EqHl limit can be reached. 

Since the volume hologram has recorded the mode 
structure of the signal beam versus reference beam inter- 
ference pattern, and since the pump beam is diffracted 
by the hologram, the diffracted beam (Ed) has the same 
mode structure than the signal beam (Et) (see EqHJ). 
This result is illustrated by Fig|3)3 where Ed is displayed 
with the same shape than Et on Fig{3^L, but with a 
smaller amplitude (77 < 1). 

The signal (Et) and diffracted (Ed) beams are thus 
spatially coherent. They can interfere constructively (or 
destructively) on a large area (~ 1 cm 2 ) light mono de- 
tector (i.e. a photodiode). This property will be useful to 
detect efficiently the tagged and untagged photons signal. 

C. Detection of the tagged photons 

The principle of tagged photons detection is illustrated 
by Fig|5l The phase up of the US beam is reversed pe- 
riodically (uo = or 7r) with period T (see FigJSt). The 
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phase of the tagged photons field E T , which follows the 
US phase, is then reversed too. 

To simplify the discussion, we will neglect the decorre- 
lation of the tagged photons field (Eqj5]is fulfilled). We 
will also modulate the phase rapidly (with respect to Tpr 
i.e. with T <C Tpr), keeping cp zero most of the time (see 
Fig|5]c), so that the hologram can be considered as static 
and unperturbed by the phase variation. 

In that case, the diffracted field Ep> will remain nearly 
constant: Ep> tj(Et)- When up is zero, Et and Ed 
are in phase, they interfere constructively and the total 
intensity signal \E T + Ep>\ 2 is maximum (see FigJSJi): 



Io = (\Et 



■ Ed\ 2 )ip=o 



\E T \ 2 (l + vf 



(6) 



When phase is tt contrarily, the Et and Ep> are opposite 
in phase, and the total intensity signal is minimum (see 
FigEb): 

J w = (\Et + E D \ 2 )^ ~ |£ T | 2 (1 - rj) 2 (7) 

Reversing the phase of the US yields to a modulation 
of the total intensity signal equal to : 



D. Detection of the untagged photons 



(8) 



It is a little bit more difficult to illustrate the detection 
of the untagged photons by a simple figure, because the 
calculation involves to consider both the untagged pho- 
tons field at the carrier frequency ujl, and the tagged 
photons fields Et and Et> which evolves at the two side- 
band frequencies ujl + for Et, and ujl — uja for Et> • 
To detect the untagged photons, we tune the reference 
beam frequency our at the untagged photons frequency: 
ur = ojl, and we modulate the US beam intensity by 
turning on and off the US beam. 

To simplify the discussion, we will neglect again the 
decorrelation of the tagged photons field (Eq|5] is ful- 
filled). We will also modulate the US beam rapidly (with 
respect to Tpr i.e. with a period T <C Tpr). Let us 
call Ejj and Ejji the untagged photons fields without, 
and with US beam. Et and Et' are the tagged photons 
fields with US (theses fields are zero without US). 

Since the energy is conserved, the total number of pho- 
tons (carrier + sidebands) does not depend on the US. 
We get thus: 



\E V \ 2 = \Eu, 



\E T \ 2 



\E T ,\ 2 



(9) 



The untagged photons field in presence of the US, e.g 
Eu>, is spatially coherent with the one without US, e.g 
Ejj- According to (Eq|9|), its magnitude can be expressed 
as follows : 



\EtjA = \EtjU 1 



\Et\ 2 + \E T >\ 2 



\Eu? 



(10) 



In practical situation, the efficiency of the acousto op- 
tic effect is low and the energy within the sideband is low 



(< 1%) with respect to the carrier. This means that the 
untagged photons field variation is low: Ejj — Ejj> <C Ejj. 
Whatever the value of the cyclic ratio modulation is, one 
can thus consider that the PR effect involves Ejj only. 
We get: 



E D ~ n Ejj 



(11) 



When the US is off, the field on the detector is Ejj + Ed 
and the detected intensity signal / is: 



\E V 



■E D \ 2 



\Ed\ 2 - 



2 V \Eu\ 2 



(12) 
(13) 



When the US is on, the field on the detector is Ejji + Ed 
for the carrier, and Et and Et> for the two sidebands. 
The intensity signal I' is: 



\Err,\ 2 - 



\Ed\ 2 



E D \ 2 + \E T \ 2 + \E T , 
- v(Eu>Eu + EuEy,) 



(14) 

\E T \ 2 + \E T ,flh) 



Taking into account the energy conservation (Eqj9j), the 
spatial coherence of (Eu,Eu>) and (EqITD|h we get the 
modulation of the detected intensity: 



i-r = 2i 1 {\e u \ 2 

~ 2tj 



\E W \.\E L 
Et\ 2 + \Et> \ 2 



2v \E T \ 2 



(16) 
(17) 



since the weight of the two sidebands components are 
approximately the same: \Et\ 2 — \Et'\ 2 • 

By comparing EqJH]and EqUHJ the detected signal have 
the same order of magnitude when detecting either the 
tagged or the untagged photons, when we consider the 
same acoustical energy. 

E. Detecting tagged or untagged photons ? 

To our knowledge, three groups are working on acousto 
optic imaging with PR detection of the signal. Two of 
them, the R.A. Royi^^-^^ and the L.V. Wang group 24 
detect the untagged photons. We are the third group 16 i 25 
and we detect both the tagged and untagged photons. 

Detection of the untagged photons is simpler since it 
is not necessary to shift frequency the reference beam 
(ujr = ujl)- The acousto optic modulators of Fig JH are 
thus not needed. Moreover, it is not necessary to ap- 
ply the US beam all the time. Untagged photons de- 
tection is thus well suited to detect very short burst of 
US beam able to give information resolved along the US 
beam propagation direction^. But short US burst yield 
a small signal, and signal is needed to image thick breast 
in vivo. 

The detection of the untagged photons corresponds to 
a small change on a large signal (white background detec- 
tion), while the detection of the tagged photons, which 
corresponds to roughly the same absolute value change, 
yields contrarily to about 100% change on a small signal 
(black background detection). Tagged photon detection 
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Fig. 6. Typical experimental setup. L: laser; BS: Beam 
Splitter; AOM1, AOM2: acousto optic modulators; PZT: 
US generator; LI, L2: light collecting lenses; PR: pho- 
torefractive crystal; PD: photo detector. 



is thus expected to give less technical noise. For example, 
vibrations on the reference beam mirrors, which modify 
the length of the pump beam arm, is expected to yield 
about 100 x 2 more noise for the untagged configuration 
than for the tagged one. 

The tagged photon configuration offers more degrees 
of freedom for the detection configuration, because the 
signal and reference beam can be modulated whether in 
phase or amplitude. 

Since we do not have tested all the possible detection 
configurations, making a complete comparison of tagged 
and untagged photons detection schemes is out of the 
scope of the present paper. We can simply say that for 
the configurations we have presently tested, the Signal to 
Noise Ratio (SNR) is about the same in the two cases. 
Since our purpose is to image breast, we need to improve 
the detection sensitivity. We continue thus to work with 
our setup that is able to detect both tagged and untagged 
photons, exploring configurations that are expected to 
yield better SNR. This work is under progress. 

4. Experimental test 

A. Setup 

A typical experimental setup, able to select either the 
tagged or the untagged photons, is shown on FigJBl The 
main laser beam is splitted in an illumination and ref- 
erence beam by the beam splitter BS. The US beam 
(2MHz) produced by the generator PZT is focused 
within the sample. The frequency offset of the illu- 
mination beam is adjusted by using the two acousto- 
optic modulators AOM1 and AOM2 which are excited 
at 7SMHz and 76MHz (for selecting the tagged pho- 
tons), or both ISM Hz (for selecting the untagged pho- 



2 here 100 is the untagged versus tagged photons field ratio 



US phase 



b) ' 




ft 
\ 














J 


V ao j 










< 1 




r 1 




■ 

4 

•f 



f . 

■ * tllTlP 



5 10 (ms) 



Fig. 7. Typical tagged photons experimental signal, a) 
Modulation of the US phase b) Tagged photons signal. 
Va.o. : maximum variation of the tagged photons signal. 
T: phase modulation period. Vertical axis is voltage at 
the output of the PD amplifier: lOOmV per div. 



tons). The light diffused by the sample is collected by the 
high NA ~ 1 (Numerical Aperture) lenses LI and L2. 
LI collects the light within the PR crystal that records 
the hologram of the selected signal beam (tagged or un- 
tagged). L2 collects the interference pattern of the signal 
beam (Et) with the diffracted beam (Ed) into the pho- 
todetector PD. 

In our setup, L is a Nd:YAG laser (1.06 /im, 1 to 5 
W CW power), the PR crystal is a 1.4 x 1.4 x 2cm 3 
GaAs crystal^, and PD is a large area photodiode (0.1 
to 0.5 cm 2 ) whose signal is amplified by a transimpedance 
amplifier (R = 100KQ to 10 Mfi). 

In the Murray's setup 17 , L is a frequency doubled 
Nd:YAG (532 nm, 80 mW), the PR crystal is a 5 x 5 x 
7mm 3 Bii2Si02o crystal, whose PR efficiency is improved 
by applying a DC electric field, the US frequency is 1.1 
MHz, and PD is an avalanche photodiode. Since Murray 
detects the untagged photons, the acousto-optic modu- 
lators are not present, but it should be pointed out that 
absorption at 532nm is more important than at 1064nm, 
and thus it can reduce the thickness of investigation. 

Our setup, which can detect both the tagged and un- 
tagged photons, is expected to be more sensitive, while 
the Murray's setup, which is used with short US pulses, 
is faster. 

B. Experimental result 

FigEl shows a typical tagged photons experimental signal 
obtained with to tt phase modulation of the US beam. 
The modulation frequency is 300Hz (modulation period 
T = 3.33 ms). The US beam frequency is 2MHz, with a 
maximum US pressure of 2MPa at the US beam waist. 
The main laser power is 1.2W. The reference and illu- 
mination beam power are both 300mVK, their areas on 
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Ultrasound position (mm) 

Fig. 8. Acousto-optic profile (a) with an optical absorber 
embedded within an Agar + Intralipid phantom of thick- 
ness t = 30mm with a reduced scatttering coefficient 
fi s = 10cm -1 @1064nm obtained with an anisotropic 
(PR) detection configuration. Light input on the sample 
is lW/cm (gaussian illumination FWHM = 1.1mm). 
The US pressure (2.3MHz,1.5MPa) is phase modulated 
at 3kHz with a duty cycle of 24%. Curve (b) represents 
a fit of the one dimensional extend of light scattering 
(e.g exp(—/jL e ffr)/r with /i e // = 2.2cm -1 ) within the US 
plane (15mm from input window), weighted by the in- 
put illumination /i e // = 2.2cm -1 . Inset represents the 
acousto-optic contrast [1 — (a) /(b)]. 



crystal and sample are both 1cm 2 . 

Measurement is performed with a 4cm chicken sample, 
whose optical properties (diffusion and absorption) are 
close to human breast. As seen, the tagged photons signal 
SNR is good (16 times averaged). One must notice that 
the signal is not rectangular like the phase modulation. 
In particular, the maximum of signal, which occurs on the 
phase plateaus (ip = tt and 0) is not flat, but decreases 
exponentially (see grey arrows on Fig|7j). 

This is expected when the PR time Tpr becomes 
shorter than the phase modulation period T. We have 
measured (see further) Tpr and we have found Tpr = 0.5 
ms. Note that Tpr can also be measured on FigJTl since 
Tpr is time constant of the grey arrow decay. 

This result is very encouraging, because it means that 
the detection bandwidth is 1/(2ttTpr) = 0.3kHz. Re- 
member the signal bandwidth is Aco = 1.5kHz (HWHM) 
on the breast. The detection is thus optimal, within a fac- 
tor 5. Since the SNR is very high (much larger than 5) 
in FigJTl we expect to get enough SNR to get significant 
result with a thick living sample. 

As an other illustration, FigJHJa) represents a pro- 
file of a Agar plus Intralipid phantom with a thick- 
ness t = 30mm and a reduced scattering coefficient 
/jl s = 1/Z* = 10cm -1 at 1064nm. The absorption co- 



efficient of water @1064nm is fi a = 0.144cm -127 . The 
sample contains an optical absorber (black ink), which 
is a cylinder with a diameter of 3mm and a height of 
7mm along the laser input direction (perpendicular to 
the US beam). The magnitude of the US pressure is 
approximately of l.bMPa at 2.3MHz, with a phase 
modulation at 3kHz and a duty cycle of 24%, corre- 
sponding to the maximum of the signal with a lock-in 
detectio n 16 ! 25 . The photorefractive holographic setup is 
based on an anisotropic diffraction configuration 28 : the 
reference beam (e.g vertically polarized) diffracts a con- 
tribution (e.g tagged-photons field) which is perpendicu- 
larly polarized (e.g horizontal); the output speckle from 
the sample is 45°-polarized from vertical direction using 
a large aperture infrared dichroic polarizer, and a sim- 
ilar analyzer is positioned in front of the photodetector 
with an horizontal polarization axis. Consequently, the 
speckle and the reference fields still interfere within the 
PR crystal in order to build the hologram, the diffracted 
reference and the speckle field recombine onto the ana- 
lyzer as well. This configuration minimizes the collection 
of the unwanted scattered reference by the PR crystal 
faces. In this experiment the tagged light is about x 10 4 
lower than the total scattered light (untagged photons 
plus scattered reference light). 

Classically, in the 3D diffusion regime and in presence 
of absorption, the spatial distribution of energy emitted 
from a point source at distance r is given by £ e -/ie//r , 

where /i e // = s/S/^ad^a + aO- This effective parameter 
indicates that attenuation is increased by scattering, that 
lengthens optical pathes. 

The continuous envelope FiglHJb) represents the fit the 
experimental data FiglHJa) using this model and taking 
into account the gaussian input illumination (FWHM = 
1.1mm). The effective coefficient /i e // is founded to 
be 2.2cm -1 , close to the theoretical value (pt e ff = 
2.1cm -1 ) given by the reduced scattering coefficient of 
the medium and the absorption coefficient of pure wa- 
ter at 1064nm defined above. The measured background 
(around 0.8mV) corresponds to the noise of the tran- 
simpedance stage of the detection, that is shot-noise lim- 
ited at this level of the scattered light. The absorbing 
element is revealed by the acousto-optic contrast, e.g 
[l-(a)/(6)], which is close to 0.22, and exhibits a FWHM 
of 7mm. This value is connected to the diameter of the 
absorbing element (3mm), the US resolution (just above 
1.5mm) and the light transport mean free path /J of the 
scattering medium (about 1mm). 

5. Measurement of the photorefractive time Tpr 

Most published results on ultrasound light modulation 
imaging have been obtained with phantoms, which do 
not decorrelate in time. In that case, the PR detection 
SNR does not strongly depends on the reference beam 
power. The power must be large enough to reach the 
plateau value for the photorefractive efficiency 77, but re- 
mains low enough to avoid noise (the reference beam is 
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Fig. 9. PR effect when the selected photons and reference 
beam frequencies are different: uur ^ uut 



and 



E D (t) oc 



Er ei Aujt f™E T (t - t) E* r e~i AujT e~ T l T ^ dr 



TpR 



\E 7 



\Er\ 



(20) 

Note that if Auu is zero, Eq|20]is identical to Eqj3l Note 
also that E R does not depend on time, and can thus be 
removed from the integral. 

Consider that a PR experiment is made with a sam- 
ple whose decorrelation time is much longer than Tpr. 
This means that decorrelation can be neglected, and that 
ET(t) is uniquely driven by the US amplitude or phase 
modulation. Erit) is thus known. \Et\ 2 can also be ne- 
glected in Eq l20l denominator (since \E T \ 2 <C l^l 2 ) and 
thus : 



scattered by the PR crystal defects yielding a parasitic 
photodiode current that brings noise). With phantoms, 
the best sensitivity is then obtained with a quite low 
power reference beam (< 100 mW in our experiment). 

With breast, the light signal is Doppler broadened 
by the tissues inner motions (brownian motion, blood 
flow...) yielding typically a spectral width of 3kHz 4 . In 
order to optimize the detection efficiency, one must in- 
crease the detection bandwidth 1/Tpr by increasing the 
reference beam power to obtain 1/Tpr ~ 3 kHz. Op- 
timal detection conditions for phantoms and breast are 
thus very different. 

To improve the detection sensitivity for future breast 
experiment, it is very important to measure Tpr. To 
get reliable result, we have proposed a technique able to 
measure Tpr in situ, i.e. in the setup that is used for 
imaging phantoms (and breast in future) 29 . 

A. Principle of the measurement of Tpr 

To perform acousto-optic imaging with living sample, it 
is very important to measure Tpr, since it is necessary 
to match Tpr with the sample decorrelation time in or- 
der to optimize the detection efficiency. The ability in 
our setup to freely choose the frequency of the reference 
beam gives new opportunity to measure the photorefrac- 
tive time Tpr. 

The idea of the measurement technique is illustrated by 
FigEl The frequency of the reference beam ujr is shifted 
with respect to the signal beam frequency uut — 
or uu L : 

our = out + Auu (18) 

To simplify the notation, we will consider that Er still 
represents a complex field projection at frequency out, so 
that Er must be replaced in the calculation by Er e^ Aujt . 
We get thus: 

{E T E* R e^ Au)t ) T 



E D (t)cx— / E T (t-T) e- jAuT e~ T/TpR dr e jAu;t 

Tpr J 

(21) 

We must notice that in Eqj2lJ Er(t) is convolved by 
two time kernels. The first kernel e~ 7 ~/ TpR is unknown 
(since Tpr is unknown), while the second e~^ AujT is 
known. Its width can be freely adjusted by tuning Aou 
with the acousto-optic modulator. 

From Eqj20] it is then straightforward to calculate the 
dependance of the acousto-optic signal with Aou for the 
different detection configurations (phase modulation for 
the tagged photons and amplitude modulation for the un- 
tagged photons)^. Comparing the calculated spectrum 
to the experiment yields then an accurate measurement 
of Tp R . 

B. Calculation of the tagged photons signal with US am- 
plitude modulation 

It is possible to calculate the tagged photons signal as 
a function of Auu in the phase modulation configuration. 
Nevertheless, as shown in^ 9 -, the shape of the spectrum 
is quite cumbersome, and it seems quite heavy to fit the 
experimental data with such a spectrum shape. 

It is thus more efficient to measure Tpr with a rectan- 
gular amplitude modulation of the US with 50% cycling 
ratio, the tagged (or untagged) acousto-optic signal be- 
ing measured with a lock-in amplifier tuned at the mod- 
ulation frequency. This is the key point of the detuning 
method, since measurements are performed in situ at the 
US modulation frequency (here 2.5kHz) and thus do not 
depend on the frequency response of the detector, which 
is quite distorted due to the many stages of electronic fil- 
ters connected to the photodetector. A straightforward 
calculation gives only three contributions for the P and 
Q quadrature of the lock-in signal 29 : 

P = P + P+ + P- (22) 

Q = Qo + Q+ + Q- (23) 

with 

24 

= i + uww <24) 



9 




-6.0 -4.0 -2.0 0.0 2.0 4.0 6.0 
Ato {kHz) 

Fig. 10. Calculated spectrum of the tagged photons sig- 
nal (amplitude modulation with 50% cycling ratio) with 
Tpr = 0.25ms (a), 0.5ms, (b), 1ms (c), 2ms (d), 4ms 
(e). Horizontal axis is the frequency offset Auj. Vertical 
axis is ^P 2 + Q 2 lock -in signal in arbitrary normalized 
units. 




-6,0 -4,0 -2.0 0.0 2.0 4,0 6,0 
Aco (kHz) 



Fig. 11. Spectrum of the tagged photons signal with 
50% duty cycle amplitude modulation of the US. Black 
curve is experimental data with reference beam flux of 
300mW/cm 2 and a modulation frequency of 2.5kHz. 
Heavy grey line curve is calculated spectrum with Tpr = 
0.45ms. Horizontal axis is the frequency offset Auj. Ver- 
tical axis is \/P 2 + Q 2 lock-in signal in arbitrary units. 



P±(Aw) 



Q±(Aw) = 



1 + {uJrnod T Aw) 2 Tj 



2 

PR 



A(w mod =F Alo)T pr 
1 + (comod T Aw) 2 T 2 R 



(25) 



(26) 



where A is a proportional constant. By using Eq|24l to 
Eql26l it is then quite simple to measure Tpp by fitting 
the experimental data with the calculated Auo spectrum. 



C. Measurement with the 
plitude modulation 



tagged photons and US am- 



From Eq|22]to Eq|26lwe have calculated the tagged pho- 
tons signal as a function of Acj, when the US beam is 
modulated with a rectangular [0, 1] amplitude modula- 
tion of 50% cycling ratio, the detection being performed 
with a lock amplifier tuned at the modulation frequency 
(2.5kHz). As seen on Fig(10]the shape of the spectrum 
is strongly dependent on Tpr. 

It is then possible to fit experimental data on the 
theoretical curves. FigfTTl shows the magnitude (e.g 
R = \/P 2 + Q 2 ) of the lock-in signal (points). The 
tagged photons are selected (cor = uol + Aco) and the 
US beam is modulated in amplitude with 50% duty cy- 
cle. The reference beam flux of 300miy/ cm 2 , and the 
modulation frequency is 2.5kHz. We have fit the ex- 
perimental data with the theoretical curve deduced from 
Eql24l to Eq|26l The fit free parameters are Tpr and 
A. The best fit yields Tpr = 0.45ms. The experimental 
data are on FigQj] as a black curve, the fit as an heavy 
grey line curve. 




0.01 



0,1 Reference Beam Flux 
(W.cm- 2 ) 



Fig. 12. Tpr in ms as a function of the reference beam 
flux in W/cm 2 . Crosses are experimental data. Dashed 
line is data extrapolation with slope —1. Grey arrow 
corresponds to FigJTl experimental conditions. 



We have recorded many spectra for different reference 
beam flux. Each spectrum has been fitted by the theo- 
retical curve yielding Tpr. FigfT2l shows in log- log scale 
Tpr as a function of the reference beam flux. Experimen- 
tal points are crosses, data linear log- log extrapolation is 
dashed line. As seen, the slope of the extrapolation line 
is —1. This means that Tpr is inversely proportional to 
the beam flux, as expected. 

The shortest photorefractive time we get is Tpr = 
0.25ms for a flux of 0.55 W/cm 2 . The FigG] modula- 
tion phase signal is obtained with Tpr = 0.5ms and 
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0.3W/ cm 2 (grey arrow on FigJTj). 
6. Conclusion 

Seeing through highly scattering media such as living tis- 
sues is a goal difficult to reach. Coupling light and ultra- 
sound in acousto-optic imaging is a promising method to 
reach this aim. Nevertheless the efficient detection of the 
tagged photons remains a challenge. 

The PR crystal detection scheme proposed here is a 
possible way to solve this problem. PR crystal detec- 
tion has many advantages. The detection optical etendue 
is large since photodetector area may be quite large 
(~ 1cm 2 ), and since the collecting lens numerical aper- 
ture can be large, too (NA ~ 1). Since the detector is 
a single-detector (photodiode), the analysis of the data 
is simple and fast. By adjusting the power of the pump 
beam, it is possible to match the detection bandwidth 
1/Tpr with the signal bandwidth Aco in order to detect 
with optimal efficiency the "tagged" or "untagged" sig- 
nal diffused by living tissues that are broadened by the 
diffuser inner motion (brownian motion, blood flow ...). 
We demonstrate here our ability to get a high SNR (see 
FigEj) with a thick chicken sample. Our chicken sample 
does not decor relate as do living tissues, the measurement 
is done with a short photorefractive time Tpr = 0.5 ms. 
This result is very encouraging. 

The results presented in this paper have been obtained 
with a Nd:YAG laser at 1064 nm and a GaAs photore- 
fractive crystal. The method could be significantly en- 
hanced by the use of a laser source at 800 nm, accord- 
ing to the absorption coefficients of hemoglobin and de- 
oxy hemoglobin, in order to perform a measurement of 
the local blood activity (two wavelengths measurements). 
We are searching at present for new PR crystals that are 
sensitive in this spectral range. 

This work is currently supported by a grant from the 
project Cancerople Ile-de-France. 
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